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Abstract 

We consider the effect of nucleation on a one-dimensional stepped surface, finding 
that step-flow growth is metastable for any strength of the additional step-edge 
barrier. The surface is made unstable by the formation of a critical nucleus, whose 
lateral size is related to the destabilization process on a high-symmetry surface. 
Arguments based on a critical nucleus of height two, which suggest the existence of 
a fully stable regime for small barrier, fail to describe this phenomenon. 
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1 Introduction 



Stability, instability and metastability are well known and widespread concepts 
in statistical mechanics, whose precise meaning may depend on the field of ap- 
plication. In this article we will refer to a crystal surface growing by Molecular 
Beam Epitaxy, which has two main growth modes: layer-by-layer growth for a 
singular, high symmetry surface, and step-flow for a vicinal, stepped surface. 
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Layer-by- layer growth [1] proceeds through diffusion of newly deposited atoms, 
the nucleation of stable islands and their increase in size via aggregation; coa- 
lescence of islands and the descent of atoms from upper terraces lead to layer 
completion. Conversely, step-flow growth [2] proceeds through attachment of 
deposited atoms to pre-existing steps. Stability means that steps keep straight 
(absence of meandering) and of uniform density (absence of step-bunching). 
The distinction between the two growth regimes is not sharp: in qualitative 
terms, step-flow dominates when the average distance between islands, the so 
called diffusion length £ D , is larger than the separation £ between pre-existing 
steps. Therefore, if m = l/£ is the slope of the surface, layer- by-layer and 
step-flow growth occur for m smaller and larger than 1/£ D , respectively. 

The presence of an additional step-edge barrier, called Ehrlich-Schwoebel 
(ES) barrier, induces a slope dependent nonequilibrium current, j E s( m ), since 
adatoms deposited on a vicinal terrace are more likely to be incorporated at the 
uphill step [3]. This barrier affects the stability properties of the surface and 
the current allows to separate in a more precise way the two growth regimes: 
from a continuum analysis [4] it turns out that a flat surface is linearly unstable 
at small slope (m < m ) and stable at large slope (m > m ), where m is the 
slope for which j ES (m) is maximal. Step-flow of a vicinal surface (m > mo) is 
then linearly stable, but this does not preclude a noise-induced destabilization 
analogous to the decay of metastable states in thermodynamics. In the present 
article we face this question: How and when a stepped surface is destabilized 
by the nucleation of mounds on vicinal terraces? 

The question, limited to the case of an infinite step-edge barrier, was faced 
by Krug and Schimschak, with Kinetic Monte Carlo (KMC) simulations [5]. 
In the limit of infinite barrier, the vicinal surface appears to be destabilized 
by the formation of islands on the terraces. This is reasonable, in particular if 
thermal detachment from steps is suppressed, but what about destabilization 
if we decrease the barrier and change the slope? 

An answer to such a question can be provided by an argument [6], based 
on the assumption that a mound of height two (i.e. a mound formed by a 
second-layer nucleation on an island) always leads to the destabilization of the 
surface. In order to understand whether the surface gets destabilized or not, 
one compares the typical time t 2nd required for second-layer nucleation, with 
the time needed for the first-layer island to be reabsorbed by the advancing 
step of the vicinal surface, t adv . Step-flow is assumed to be destabilized only 
if t 2nd < t adv . Explicit evaluation of these quantities indicates that the time 
t 2nd increases as the ES barrier goes to zero, while t adv is reduced. It turns 
out then [7] that the criterion implies the existence of a minimal barrier for 
destabilization and, as a consequence, that step-flow growth should be fully 
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stable at small barriers. 1 

In this Letter we show that the above picture is not correct. First, we present 
numerical evidence that the time after which the system becomes unstable 
is finite for all finite values of £ ES , and it diverges (exponentially) only for 
£ ES — ► 0. Then, we provide and validate a heuristic argument that elucidates 
the destabilization mechanism. 



2 The destabilization process 

We perform Kinetic Monte Carlo simulations in one dimension of the simplest 
model taking into account deposition, diffusion and the presence of an inter- 
layer ES barrier. The surface is represented by a set of integer height variables 
hi on a lattice of L sites. An average tilt m is imposed via helical boundary 
conditions. The initial condition is a regular train of steps, separated by a 
terrace size £ = 1/m. Deposition events occur at rate F on randomly selected 
sites. Adatoms attempt diffusion hops to nearest neighbor sites at rate D if 
the neighbor belongs to the same terrace, and at rate D' < D if the atom 
must descend a step. Dimers and larger islands are immobile, and no thermal 
detachment of atoms from steps is allowed. We take F — 1 and D = 5x 10 5 , 
so that the diffusion length is £ D 40 (the same value used in [5]). The ES 
length can be written as £ ES = D/D' — 1 [4]. 

In order to characterize quantitatively the destabilization process we focus 
on the formation of the first "critical nucleus". This concept arises in the 
continuum theory of metastability as a localized stationary solution of the dy- 
namics, separating stable from unstable solutions. In order to explain how we 
have identified it numerically, let us briefly sketch how destabilization comes 
about (Fig. 1). 

During the early stages of growth the surface is made of a train of vicinal 
terraces with lengths fluctuating around the average value £ = 1/m. Once in a 
while new islands are nucleated on them. Some of them may grow into small 
bumps that are later reabsorbed by advancing steps. Destabilization occurs 
when a bump reaches some critical size (critical nucleus) and it starts to grow 
irreversibly. In principle the destabilization time r is the moment when mound 
growth becomes irreversible. However, in a simulation such a moment is not 



1 If we introduce the ES length £ E3 = exp(AE/kBT) — l, where AE is the additional 
step-edge barrier and T is the temperature, the minimal barrier corresponds to the 
critical length £ ES ~ cq£^/£ 3 , with a prefactor Co > 10. Taking £ D ~ 40 (valid for 
our simulations, see below), we get, for m = 0.1, £ ES > 10 4 , while we numerically 
see metastability down to £ ES = 12. 
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Fig. 1. Schematic plot of the growing surface at different times, illustrating a 
typical profile during the intermediate metastable regime (bottom), followed by the 
nucleation of a dimer and the formation of a supercritical mound. 

well defined, because mound growth is a stochastic process and even mounds 
which are not reabsorbed may temporarily shrink to some extent during their 
growth, making the numerical distinction between stable and unstable mounds 
very hard. 

We decided to set a height threshold at 20 layers: when a mound reaches such 
a height we consider step-flow to be destabilized and measure the destabiliza- 
tion time r. There are two possible systematic errors in this way of proceeding: 
in the first place, one may take as unstable a mound that is going to be re- 
absorbed; secondly, one overestimates r, since the true destabilization of a 
mound typically happens before it gets 20 monolayers high. We have carefully 
taken into account these possible problems. For all values of £ ES and m con- 
sidered, we have never seen mounds higher than 14 layers being reabsorbed. 
It is therefore very likely that all mounds we have taken as unstable (higher 
than 20 layers) were actually so. For what concerns the determination of r, 
the time needed for a mound to grow from height 1 to 20 is of the order of 
few tens of monolayers. This may introduce a significant systematic error and 
invalidate the results when r is small, but it is irrelevant when destabilization 
occurs after thousands of monolayers or more, which is the typical case here. 

The destabilization process can be split into: (1) the formation of a dimer on 
a vicinal terrace (which occurs at a rate l/r nuc j(L)); (2) the evolution of the 
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Fig. 2. Double logarithmic plot of the logarithm of l/p u = t{L) /r nuc /(L) as a 
function of the ES length for three values of the slope m, computed for a sys- 
tem of size L = 1000 and for the RW model for m = 0.2. Time is measured in 
monolayers, i.e. 1/F is the time unit. Statistical error bars are of the order of the 
symbol size. Least-square fits (shown as lines) to the functional form exp(a/iZs) 
yield 7 = 0.49,0.51 and 0.51 for m = 0.1,0.2 and 0.5, respectively. 



dimer into an unstable mound (occurring with probability p u ) . The total rate 
for the destabilization is the product of these two quantities, so that for a 
system of size L the mean destabilization time r(L) is 

r(L) = Tnud (L)-. (1) 

Pu 



r nud{L) goes to a finite value for £ ES — > and contains the whole 1/L depen- 
dence on the system size. p u instead, does not depend on L and contains the 
relevant dependence on £ ES . 

In Fig. 2 we plot, in a double logarithmic scale, the logarithm of l/p u as a 
function of the Ehrlich-Schwoebel length £ ES for several values of the tilt m. 
It is clear that the instability time diverges for £ ES — ■> as exp(a/£1 s ); there 
is no indication of a divergence at finite £ ES , i.e. no indication of a fully stable 
region. The value of 7 is of the order of 0.5 for all slopes considered. For large 
£ ES the curves approach a constant value increasing with m. 
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Fig. 3. Double logarithmic plot of k max as a function of the ES length for three 
values of the slope m. Least-square fits (shown as lines) to the functional form i^s 
yield (3 = 0.37 ± 0.01 for m = 0.1, (3 = 0.39 ± 0.02 for m = 0.2 and = 0.39 ± 0.03 
for m = 0.5. 

3 The critical nucleus 



In order to look for a theoretical explanation of the behavior of l/p u and 
understand why the argument based on a critical nucleus of height two [6] 
fails, it is important to get information on the shape and the size of the critical 
nucleus. For this purpose, we compute for each realization i of the dynamics, 
the height k max ,i of the highest mound that is eventually reabsorbed without 
diverging. After averaging over many realizations we obtain the values k max 
that are reported in Fig. 3. In the hypothesis that a critical height k c can be 
defined, k max is clearly a lower bound for k c . Its growth with decreasing £ ES 
invalidates the assumption k c = 2 upon which the argument of Ref. [6] is built: 
Figure 3 shows that, in order to become unstable, mounds must grow higher 
and higher as £ ES is reduced. 



We now relate these observations to the exponential divergence of the insta- 
bility time, via a reformulation of the mound dynamics as a one-dimensional 
random walk (RW). The formation of a mound begins with the nucleation of a 
dimer on a vicinal terrace of size £ occurring at a rate l/r nuc i(L). After a dimer 
is formed, two possible events may occur: either the dimer is reabsorbed, with 
probability 1 — p + , or it grows higher, i.e. a second- layer nucleation occurs with 
probability p + before the advancing step catches the dimer. In the latter case 
a mound of height 2 is formed and the ensuing evolution may be the decay 
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into a mound one-layer high or the upgrade to a mound of height 3. Hence the 
height A; of a mound performs a one-dimensional random walk. Particles are 
introduced at k — 1 via dimer nucleation on a terrace and may disappear at 
a sink in k = 0. The boundary condition for high k encodes the way a mound 
becomes unstable. It is possible to extract from the KMC simulations the de- 
tailed form of p+(k) [7] and it is found that p+(k) grows with k from to a 
constant value slightly larger than 1/2 for k > k s , with k s growing as £ ES — > 0. 
Hence a critical height (k s ) can be defined only in a probabilistic sense: the 
mound becomes unstable because once it reaches, through a rare fluctuation, 
the height k s it is on average more likely for it to grow higher than to shrink. 
Using the precise form of p+(k) it is possible to determine numerically \jp u 
within the RW model. The result, presented in Fig. 2 shows an excellent agree- 
ment, indicating that the RW reformulation accurately describes the mound 
evolution. 



As just discussed, it is not fully appropriate to consider a deterministic critical 
height k c . However, if we approximate the true form of p+{k) with a value 
p° + < 1/2 for k < k c and with 1 for k > k c the probability p u can be easily 
determined [8] 



IT 
i-p° 

~2T 



(2) 



Since p° + is smaller than 1/2 and k max (related to the critical height by k max = 
k c — 1) is seen to increase in the limit £ ES — > (Fig. 3), then 

l/Pu ~ (1/P+ - l) fcc ~ exp(a'k max ), (3) 



with a' depending on p+, i.e. on the slope m. Despite the approximations used 
in its derivation, formula (3) gives a divergence exp(l/^? s ) of the instability 
time, with (3 of the order of 0.4 (see Fig. 3), in rather good agreement with 
the KMC simulations. 

In order to get deeper insight into the destabilization mechanism, let us con- 
sider the width A c of the critical nucleus. To measure it in our simulations an 
unambiguous operational definition of the width A of a mound is needed. The 
simplest would be to take A as the base of the mound, but nucleation events 
occurring along the sides make difficult a precise identification of the base. 
As a consequence of that, reasonable definitions of the base lead to strongly 
fluctuating quantities. We have found that a convenient way to proceed is to 
take A as the size of the sub-top terrace, i.e. the second highest terrace in 
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Fig. 4. Double-logarithmic plot of the critical mound width A c as a function of the 
ES length for three values of the slope m. Statistical error bars on A c are of the 
order of the symbol size. As a guide for the eye, a divergence £ E ^ 2 is also plotted. 

the mound. We have considered critical mounds 2 and determined the corre- 
sponding critical width A c , as a function of £ ES for three different values of to 
(Fig. 4). 

The most important feature of Fig. 4 is that there is a unique curve A C (£ ES ) for 
all values of to: when the mound becomes unstable, the width of the sub-top 
terrace depends only on £ ES and not on the average surface tilt. Moreover, we 
find that A C (£ ES ) grows with decreasing £ ES as a power-law with exponent not 
far from 1/2. 

This phenomenology can be consistently interpreted by assuming that a mound 
formed on a vicinal surface is similar to a mound appearing on a singular sur- 
face of size L of the order of A. Let us recall that for a high symmetry surface 
there exists a critical size L C (£ ES ) such that a surface larger than L c is lin- 
early unstable, while, if the opposite is true, the instability is hindered and 
the planar surface is stable [9]. For mounds on stepped surfaces this threshold 
is reflected in the existence of a critical width A C (£ ES ), close to, although not 
necessarily coinciding with L C (£ ES ). As long as a mound on a stepped surface 
is small it tends on average to be reabsorbed, exactly as perturbations on a 
high-symmetry surface with L < L c tend to disappear. However, if nucleation 
events occur more rapidly than on average, so that its lateral size becomes 



2 We have considered the mounds that reached the threshold height 20 and taken 
as A c the size of their sub-top terrace when their height was 10. 
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larger than 

A-c(^es) ~ L C (£ ES ), 



(4) 



then the mound becomes unstable. 

This scenario is fully consistent with the numerical results presented in Fig. 4: 
in particular, the independence of A C (£ ES ) on m and the divergence of A C (£ ES ) 
for £ ES — > (consistent with the fact that, in a continuum treatment, L c m 
1/v^es) strongly support the validity of Eq. (4). 

The connection between the instability on a high-symmetry surface (which 
occurs independently of the amplitude of fluctuations) and on a stepped sur- 
face has an additional implication: the relevant quantity is the width and not 
the height of the mound. The observed increase of the critical nucleus height 
as £ ES — > 0, is a consequence of the more fundamental growth of its lateral 
size. We have numerically tested this by taking mounds larger than A C (£ ES ) 
and drastically reducing their height to 2 or 3 layers. This modification does 
not change their unstable nature. This observation confirms that the height 
of a mound is not the relevant variable for the instability. Nevertheless, the 
width and the height of a mound are dynamically correlated and for this rea- 
son studying the height of a mound is useful even if its stable or unstable 
character depends on its width. 



4 Conclusions and open questions 

In summary, we have shown that any additional step-edge barrier makes step- 
flow growth unstable with respect to the formation of mounds. The time 
needed for this instability to take place diverges exponentially as £ ES goes 
to zero. This is a consequence of the divergence of the size of the critical nu- 
cleus needed to destroy the metastable step-flow. This divergence is, in its 
turn, related to the process of destabilization of a high-symmetry surface. 

There are several open, or not fully understood, questions which we are cur- 
rently investigating [7]. The first issue is to consider how our arguments can be 
extended to two dimensions. When steps are lines rather than points, the addi- 
tional possibility of meandering arises and, as shown by Kallunki and Krug [6], 
this is a second mechanism for destabilizing step-flow, which is effective even 
if nucleation is completely suppressed. As for the destabilization induced by 
nucleation, the k c = 2 approximation still provides [6] a threshold £ ES . Since 
this result is based on a constant critical height k c , independent from £ ES , it 
is most likely that the approximation fails in two dimensions as well. 
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With regard to numerics, we want to analyze the crossover between the lin- 
ear unstable regime (m < m ) and the metastable regime (m > m ). This 
crossover is made subtle by the £ ES — dependence of the slope m . In addition, 
we want to study the effect of thermal detachment from steps. Does it have 
purely quantitative effects or does it change in a profound way the destabi- 
lization picture? Finally, from the analytical point of view, we have found that 
standard continuum (Cahn-Hilliard type) theories are not suitable for describ- 
ing the £ ES dependence of the destabilization process of the vicinal surface, not 
even in a qualitative way [7]. The possibility to give an alternative continuum 
description is an important open issue. 
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